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Using gamma-ray data from observations of the Milky Way, Andromeda (M31), and the cosmic 
background, we calculate conservative upper limits on the dark matter self-annihilation cross section 
to monoenergetic gamma rays, (aAv)-yy, over a wide range of dark matter masses. (In fact, over 
most of this range, our results are unchanged if one considers just the branching ratio to gamma 
rays with energies within a factor of a few of the endpoint at the dark matter mass.) If the final- 
state branching ratio to gamma rays, Br^y), were known, then (a Av)- n / Br (77) would define an 
upper limit on the total cross section; we conservatively assume Br(77) > 10 -4 . An upper limit 
on the total cross section can also be derived by considering the appearance rates of any Standard 
Model particles; in practice, this limit is defined by neutrinos, which are the least detectable. For 
intermediate dark matter masses, gamma-ray-based and neutrino-based upper limits on the total 
cross section are comparable, while the gamma-ray limit is stronger for small masses and the neutrino 
limit is stronger for large masses. We comment on how these results depend on the assumptions 
about astrophysical inputs and annihilation final states, and how GLAST and other gamma-ray 
experiments can improve upon them. 

PACS numbers: 95.35.+d, 95.85.Pw, 98.70.Vc, 98.62.Gq 



I. INTRODUCTION 

While there is ample gravitational evidence for dark 
matter, the nature of these particles remains mysteri- 
ous and is defined principally by the weakness of their 
interactions; for reviews, see e.g., Refs. 0, |2, The 
dark matter self-annihilation cross section is of particular 
importance, since dark matter concentrations will pro- 
duce gamma rays and other detectable Standard Model 
(SM) particles. If the dark matter (DM) is a thermal 
relic of the early universe, the annihilation cross section 
must be (<jav) ~ 3 x 10~ 26 cm 3 s~ 4 in order to obtain 
the observed relic abundance, Odm — 0.3. (Through- 
out, we consider this cross section averaged with veloc- 
ity over the dark matter velocity distribution; in the 
Milky Way, v rms ~ 10~ 3 c.) It is possible that dark 
matter is not a thermal relic, e.g., Ref. 0, [![, which 
makes it even more interesting to consider direct late- 
universe constraints on the annihilation cross section, 

e.g. , Refs. 0, a 1 i 0, s, a m m m . 

Even if the total annihilation cross section is set by 
the relic abundance, the branching ratios to specific fi- 
nal states are model-dependent. The dark matter dis- 
appearance rate due to annihilation can be constrained 
by the appearance rates of various SM particles. If the 
dark matter is the lightest stable particle in some new 
physics sector, then it can be natural to have the final- 
state branching ratio to SM particles, Br(SM), be 100%, 
as we assume. If the final states include new and purely 
sterile particles, then all appearance-based results are 
weakened proportionally to Br (SM). 

We assume that annihilation is not prevented in prin- 
ciple by dark matter not being its own antiparticle, or in 



practice by a large particle-antiparticle asymmetry. We 
also assume that a single type of new particle comprises 
the dark matter that is required to exist in the present- 
day universe, and that, consistent with observations, the 
density distributions of dark matter halos are not ap- 
preciably affected by possible dark matter interactions. 
These assumptions are made implicitly in nearly all pa- 
pers about dark matter annihilation. 

Here we calculate the constraints that can be placed 
on the annihilation cross section using gamma rays, the 
most detectable final states, over a wide range of dark 
matter masses. We first focus on the 77 final state, as it 
would be a very clean signature of dark matter annihila- 
tion, with E 1 — m x , e.g., Refs. [HI, 

Unfortunately, in typical models, this is small, Brjj'y) ~ 
1CU 4 — 10~ 3 ; in some models, it can be larger [20, l2ll. |22| . 
but one cannot be certain that these predictions match 
with nature. Since gamma rays will be ubiquitously pro- 
duced, directly in SM final states, or through radiative 
corrections and energy-loss processes, we also consider 
more general outcomes, in which the gamma-ray ener- 
gies are in a broader range below m x . 

We consider constraints on the dark matter annihila- 
tion cross section over a large mass range of 10 -5 - 10 5 
GeV. At all but the highest energies, gamma-ray data is 
available to test the annihilation cross section, provided 
that we combine constraints defined using the Milky Way 
halo, the Andromeda halo, and all the halos in the uni- 
verse. (Modern data, especially those from observations 
of the Milky Way and Andromeda, are significantly more 
constraining than those that were available earlier, e.g., 
at the time of Ref. [23], which considered limits on the 
decay of an unstable massive neutrino.) We hope that 
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our results will be useful in challenging experiments to 
report stronger limits using new data and focused analy- 
ses. With the launch of GLAST this year, and with new 
studies by TeV-range experiments, these prospects are 
good. Using our upper limits on the dark matter anni- 
hilation cross section to gamma rays, and a conservative 
assumption about the branching ratio to monoenergetic 
gamma rays, we define upper limits on the total cross 
section and compare to other constraints. 

Since the dark matter annihilation rate scales with 
density squared and the density profiles are uncertain, 
we are mindful of how our constraints on the cross sec- 
tion are affected by astrophysical uncertainties. We are 
conservative in our input choices and analysis methods, 
and we show how our results depend on these. In light of 
these considerations, we do not consider corrections that 
would affect the results by less than a factor of ~ 2, which 
also allows some simplifications. Our upper bounds on 
the annihilation cross section to gamma rays would only 
be improved by more optimistic assumptions. 

In Section II, we discuss important general bounds on 
the total annihilation cross section. In Section III, we 
review the analysis methods used for the case of gamma- 
ray lines from various dark matter concentrations. The 
experiments and observations we use are discussed in Sec- 
tion IV. In Section V, we summarize and interpret our 
results. 



nificantly distorted by dark matter annihilation. Note 
that this limit is very weak for all but the lightest masses. 



We now discuss limits which arise from the appear- 
ance rate of dark matter annihilation products, assum- 
ing Br(SM) = 100%. All final states except neutrinos 
obviously produce gamma rays, either directly or as sec- 
ondary particles (we return to neutrinos next). Quarks 
and gluons hadronize, producing pions and thus photons 
via 7T° — > 77, while r^, IV ± and Z° also produce n° 
via their decays. Charged particles produce photons via 
electromagnetic radiative corrections [271 ] . and electrons 
and positrons also produce photons via energy loss pro- 
cesses 28]. Therefore we expect a broad spectrum of 
final state photons, even though the branching ratio to 
the monoenergetic 77 final state may be small. We use 
the gamma-ray data to place upper limits on (ctav)^, 
and these are of general interest for their own sake. With 
an assumption about Br (77), these results also define an 
upper limit on the total cross section of (aAv)'yy/Br("f"f). 
Unless one is confident that Br (77) cannot be too small, 
then this constraint on the total cross section can be ar- 
bitrarily weakened. 



II. CROSS SECTION CONSTRAINTS 

The annihilation cross section sets the dark matter 
disappearance rate, for which there are two important 
constraints. The first is unitarity 0, [25|, which sets 
a general upper bound that can only be evaded in un- 
usual cases [261 ] . In the low- velocity limit, where s-wave 
annihilation dominates, the unitarity bound is {o~av) < 
An/rri^v, or 



(a A v) < (1.5 x 10 



" 13 cm 3 s- 1 ) 



300 km/s 



GeV 



(1) 

For m x > 10 6 GeV, this would require that {<tav) be 
smaller than that for a thermal relic, and which would 
thus give too large of a relic abundance [24j |. However, 
for smaller m x the unitarity limit is much less constrain- 
ing. The second general constraint comes from the re- 
quirement that annihilation does not drastically alter the 
density profiles of dark matter halos in the Universe to- 
day. In the model of Ref. [4(, a large self-annihilation 
cross section was invoked in order to reconcile predicted 
cuspy density profiles with the flatter ones inferred from 
observation, requiring 



(ctav) 



KKT 



(3 x 10" 19 cm 3 s -1 ^ 



GeV 



(2) 



We re-interpret this result as an approximate upper 
bound, beyond which halo density profiles would be sig- 



An important general limit on the total annihilation 
cross section can be obtained by considering annihila- 
tion into the least detectable final state, namely neutri- 
nos @, [1] • Given that stronger constraints will exist on 
all final states other than neutrinos, we can set a con- 
servative upper bound on the total dark matter annihi- 
lation rate by assuming the branching ratio to neutri- 
nos is 100%. (Unlike all other constraints, the neutrino 
constraint, being the weakest, is not to be divided by a 
realistic branching ratio; this follows from the fact that 
the sum of all branching ratios must be 100%.) The 
resulting limits are surprisingly strong. Dark matter an- 
nihilation into neutrinos was explored by Beacom, Bell, 
and Mack (BBM) 0, and Yiiksel et al. (YHBA) [|, for 
cosmic and Galactic dark matter sources, respectively. 
By requiring that the neutrino flux produced by anni- 
hilation be smaller than the measured atmospheric neu- 
trino background, robust bounds on the total annihila- 
tion cross section were obtained over a wide mass range. 
(Ref. 2^] extended the dark matter annihilation limits 
to lower masses and Ref. [3(| developed analogous dark 
matter decay limits over a wide range of masses.) For all 
masses considered, these limits are much stronger than 
the KKT limit; they are also stronger than the unitar- 
ity limit except at high masses. While neutrinos are the 
least detectable annihilation products, even they are ac- 
companied by gamma rays via electroweak radiative cor- 
rections; these results lead to constraints on (o~av) that 
are comparable to or better than those obtained directly 
with neutrinos [1, [M [10] ■ 
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III. CALCULATION OF DARK MATTER 
SIGNALS 

The dark matter annihilation rate depends on the 
square of the dark matter number density p/m x , which 
is written in terms of the unknown mass m x and the 
uncertain dark matter mass density p. Not coinciden- 
tally, where the density is largest, at the centers of halos, 
the uncertainties are the largest; these regions contribute 
relatively little to the gravitationally-measured mass of 
a halo. To cover as large of an energy range as possible, 
we have to consider gamma-ray data for the Milky Way, 
Andromeda, and all of the dark matter halos in the uni- 
verse. In all cases, though the astrophysical and analysis 
uncertainties vary in their severity, we make conservative 
choices for the dark matter density and hence the cross 
section limits (smaller choices for the density mean larger 
upper limits on the cross section). 

A. Dark Matter Halos 

A standard parameterization of the dark matter den- 
sity profile in a halo is 

P[r) = (r/r s r[l + (r/r s rf~^ ' ^ 

The Navarro-Frenk- White (NFW) [H[ and Kravtsov pro- 
files [32j are defined by (a, /?, 7) = (1, 3, 1) and (a, /?, 7) = 
(2,3,0.4) respectively. Near the center of the halo, the 
density of an NFW profile scales with radius as 1/r, 
while the Kravtsov profile scales less steeply, as 1/r . 
For large radii, r > r s , these two profiles coincide more 
closely. For the Milky Way, r s is 20 kpc for the NFW 
profile and 10 kpc for the Kravtsov profile, while the 
normalization, po, is fixed such that the density at the 
solar circle distance, R sc — 8.5 kpc, is p(R sc ) = 0.3 GeV 
cm" 3 (0.37 GeV cm" 3 ) for the NFW (Kravtsov) profile. 

There does not appear to be a consensus on the val- 
ues of the halo parameters for Andromeda; for example, 
compare the NFW profiles in Ref. 33[ with Ref. 
where both po an d r s are quite different. Thus we have 
chosen to model Andromeda using the Milky Way param- 
eters, as an appropriate compromise between competing 
extremes. 

In the innermost regions of halos, the uncertainties in 
the dark matter density and thus the annihilation rate 
are at their largest. For larger central regions, or whole 
halos, these uncertainties are much less. For the Milky 
Way, these effects can easily be seen in Fig. 2 of YHBA; 
for larger angular regions centered on the Galactic Cen- 
ter, the dark matter annihilation signals for different pro- 
files become much more similar to each other. While 
the uncertainties at small angular scales can be orders 
of magnitude, those at large angular scales are not more 
than a factor of about 2. 

In YHBA, we explored in detail how various annihila- 
tion signals depend on the choice of dark matter density 



profile. Our overall approach is to be conservative by 
adopting smaller choices of the astrophysical inputs; this 
means that larger values of the cross section would be 
required to get the same gamma ray or neutrino fluxes. 
Here we use the Kravtsov profile for our main results; 
for the commonly-adopted NFW profile, we find smaller 
(more restrictive) upper bounds on the cross sections, as 
shown below. Also to be conservative, we neglect the pos- 
sibility of halo substructure, e.g., Refs. [35l l36l|. or mini- 
spikes around intermediate-mass black holes [13, [H[, 
which would lead to enhanced annihilation signals. 



B. Milky Way and Andromeda Signals 

We first consider annihilations in our Galaxy, follow- 
ing the conventions of YHBA, and generalize this to the 
nearby galaxy Andromeda (M31). The intensity (flux 
per solid angle) of the annihilation signal at an angle ip 
with respect to the Galactic Center (GC) is proportional 
to the square of the dark matter density integrated over 
the line of sight, 

J(ip) = J J p 2 (V# s 2 c -21R sc cosV> + ^ 2 ) dl , (4) 

where Jo = 1/ [8. 5 kpc x (0.3 GeV cm" 3 ) 2 ] is an arbitrary 
normalization we use to make J a dimcnsionless quan- 
tity, and which cancels in the final results. The upper 
limit of the integration is given by £ m ax = (-Rjyrw — 
sin 2 tpR^) 1 ' 2 + R sc cosip. We define Jc&i as the average 
of J over a cone of half-angle tp centered on the GC, 

Jaq = ^iJ q sin W , (5) 

where A57 = 27r(l — cos?/') is the angular size of the cone 
in steradians. The values of i7(V0 and Jmi can be read 
directly from Fig. 2 of YHBA (below, we do not explicitly 
show the sr _1 units of J{ip) and J&n)- 

Eq. ([41 can easily be generalized to external ha- 
los [39|, |4fj| (such as the Andromeda galaxy at a distance 
of I?M3i — 700 Mpc) using 

J(i>) = Jo J P 2 {^J £&31 - 2 ^M31 COS ifl + A di , 

where the result is independent of the upper and lower 
limits of integration (l m in,£max) as long as they cover 
most of the halo under consideration. 

For extragalactic dark matter sources, the annihilation 
signals will include a contribution from the dark matter 
in our own galaxy along the line of sight. However, in the 
case of an external galaxy like Andromeda, this contri- 
bution will be eliminated if there is a subtraction of the 
background intensity from a region close to the source, 
as is often done in observational analyses. 
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With these definitions, the intensity of the dark matter 
annihilation gamma-ray signal is 



d*&2 _ (<jav) Jaq 1 dN 1 
dE ~ 2 J 47rm 2 dE 



(7) 



where dN^/dE is the gamma-ray spectrum per annihila- 
tion. In the case of annihilation into two monoenergetic 
gamma rays, we simply have dN^/dE = 25(m x — E); we 
generalize this below. Similarly, the total flux (per unit 
energy) from a region of solid angle AS! is 



d(p 



2 _ ^i An _ (PAV) Ja^AS! 1 dN. 
dE ~ dE ~ 2 J 47rm 2 ^E 
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(8) 



C. Cosmic Diffuse Signal 

The calculation of the cosmic diffuse annihilation sig- 
nal is detailed, for example, in Rcfs. @, |4l|, where the 
cosmological flux of annihilation products from external 
galaxies was calculated taking the clustering of dark mat- 
ter into account. Recently, BBM and YHBA applied this 
technique to the scenario where dark matter annihilates 
into neutrinos. 

The cosmic diffuse flux, arising from dark matter an- 
nihilation in halos throughout the Universe, is 



d$ 7 
~~dE 



O 2 n 2 



m 2 



■f(z)(l + z) 3 dN,(E>) _, 



h(z) 



dE' 



w *dz, (9) 



where Hq = 70 km s~ x Mpc^ 1 is the Hubble parame- 
ter and O.DM i s ^ ne dark matter density in units of the 
critical density. We assume a flat universe, with SIdm 
= 0.3, n A = 0.7, h{z) = [(1 + zfn DM + S! A ] 1/2 . The 
factor e -z / Zmax , taken from Ref. p, accounts for the at- 
tenuation of gamma rays, a modest effect for the energies 
considered here. The factor f(z) in Eq. ^ accounts for 
the average increase in density squared due to the fact 
that dark matter is clustered into halos, rather than uni- 
formly distributed, and the evolution with redshift of the 
halo number density. (The A 2 factor in BBM is equal to 
f(z)(l + z) 3 .) Following YHBA, we use the parameteriza- 
tion log 10 (/(z)// ) = 0.9 [exp(-0.9z) - 1] - 0.16z, where 
/o depends on the halo profile. Choosing the Kravtsov 
(NFW) profile, / ~ 2 (5) x 10 4 . 

Gamma rays that are produced with energy E' are 
observed with rcdshiftcd energy E = E'/(l + z). For 
annihilation into monoenergetic gamma rays, the delta 
function source spectrum is modified by redshift as 



dN „ w , 2 . / ,m y 

= 2 5(m y - E') = — 8(z- (— - 1 

dE' y x E \ E 



(10) 



which shows that the observed flux at an energy E is 
contributed by sources at redshift — 1. 



IV. SPECIFIC OBSERVATIONS AND DERIVED 
ANNIHILATION CONSTRAINTS 

We have collected gamma-ray flux measurements and 
limits from a wide variety of experiments, spanning an 
extensive energy range from 20 keV to 10 TeV. In most 
of the observations, the energy spectra are given in 
log-spaced energy intervals. We calculate annihilation 
gamma-ray fluxes for the Galactic, Andromeda, and cos- 
mic dark matter sources, using the methods outlined in 
Section III above. These are compared with observa- 
tional data over an energy range, conveniently chosen as 
10~°' 4 m x - m x , that is comparable to or larger than the 
energy resolution and bin size of the experiments. If only 
upper limits on the flux are given, we instead compare 
our predictions directly with these upper limits. 

Our constraints on the dark matter annihilation rate 
are conservatively determined by demanding that the 
annihilation flux be smaller than 100% of the observed 
(presumably not produced by dark matter) gamma-ray 
background flux at the corresponding energy range. In 
Fig. [I] we show the GC and cosmic diffuse signals from 
dark matter annihilations which fulfill this criterion, su- 
perimposed upon the Galactic and extragalactic spectra, 
respectively, as measured by COMPTEL and EGRET. 

The experiments report their results as either intensity 
(as in Eq. [7]), which requires that we calculate J7asi, or 
flux from a given angular region (as in Eq. [5J , for which 
we need J'j^qAD,. We present the values of these param- 
eters which correspond to the Kravtsov profile, as this 
results in the most conservative upper limits on the an- 
nihilation cross section. Our limits on the dark matter 
annihilation cross section are reported in Fig. [21 where 
we also show how our results would change if the NFW 
profile were adopted instead. The details of the experi- 
ments and our analyses are summarized below for each 
observation. 



A. COMPTEL and EGRET 

COMPTEL 42], the imaging COMPton TELescope 
aboard the Compton Gamma Ray Observatory (CGRO) 
satellite, measured gamma rays in the energy range 1— 
30 MeV. EGRET [43], the Energetic Gamma Ray Ex- 
periment Telescope, also aboard the CGRO, measured 
gamma rays in the energy range 30 MeV to nearly 100 
GeV. For both COMPTEL and EGRET, the full sky was 
studied with an angular resolution of at worst a few de- 
grees (for the large regions we consider, this makes no 
difference). The energy resolution was modest, and the 
data were given in a few logarithmically-spaced bins per 
decade in energy. 

Both COMPTEL and EGRET observed the Galactic 
Center region, and the measured gamma-ray intensity 
energy spectra are reported in Ref. [44j, |45| for the region 
-30° < I < 30° and -5° < b < 5° (Galactic longitude 
and latitude, respectively). The disk-like morphology of 
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E [ GeV ] 

FIG. 1: Example dark matter annihilation signals, shown 
superimposed on the Galactic and extragalactic gamma-ray 
spectra measured by COMPTEL and EGRET. In each case, 
the cross section is chosen so that the signals are normalized 
according to our conservative detection criteria, namely, that 
the signal be 100% of the size of the background when in- 
tegrated in the energy range chosen (0.4 in log 10 E, shown 
by horizontal arrows). The narrow signal on the right is the 
Galactic Center flux due to annihilation into monoenergetic 
gamma rays, for m x = 1 GeV; the signal is smeared as appro- 
priate for a detection with finite energy resolution. The broad 
feature on the left is the cosmic diffuse signal for annihilation 
into monoenergetic gamma rays at m x =0.1 GeV, smeared 
by redshift. 



the emission region makes it clear that nearly all of this 
emission is due to ordinary astrophysical sources; to be 
conservative, we do not attempt to define a limit on the 
small component of this that could be due to centrally- 
concentrated dark matter, and simply use the total ob- 
served intensity to bound any dark matter contribution. 
Also, we evaluate the dark matter signal as if from a cir- 
cular region of ip = 30°; accounting for the rectangular 
shape of the region would lead to a higher value than the 
J An — 13 that we adopt. Using a less conservative set of 
assumptions than we employ, stronger limits on (<jv) 11 
were derived from the EGRET data in Ref . [H[ . 



B. H.E.S.S. 

H.E.S.S. (High Energy Stereoscopic System), a system 
of multiple atmospheric Cerenkov telescopes, is presently 
in operation in Namibia [4(|. H.E.S.S. observed the 
Galactic Center region in the energy range 0.3-15 TeV. 
An apparent point source at the Galactic Center was ob- 
served, as was an extended source (~ 1°) known as the 
Galactic Center Ridge [13] ■ While the origin of the point 



source is unknown, the Ridge emission is almost certainly 
astrophysical, and is consistent with being caused by cos- 
mic rays colliding with a gas cloud (again, we do not at- 
tempt to account for this, and will simply bound any dark 
matter contribution by the total observed intensity). 

Since the uncertainties in the dark matter profile in- 
crease for smaller angular regions around the Galac- 
tic Center, it is more robust to define our results us- 
ing the extended region instead of the point source. 
The Ridge emission was observed in an angular region 
-0.8° < I < 0.8° and -0.3° < b < 0.3° in Galactic 
coordinates, and the resulting flux reported by H.E.S.S. 
reflected a background subtraction from a nearby region 
(-0.8° < I < 0.8° and 0.8° < b < 1.5°) to help ac- 
count for cosmic rays. Thus, we have to consider not the 
whole dark matter signal, but just its contrast between 
the central and adjacent regions by accounting for this 
subtraction in our analysis. 

We approximate the intensity from the rectangular re- 
gion of the Galactic Center Ridge with a circle of radius 
0.8°. We also estimate the adverse effect of the back- 
ground subtraction on our limits by choosing J to be 
subtracted at its maximum, i.e., ip = 0.8°. This means 

n rO.S" 

jAQ = MlJ (W)-^(°- 8 °)) sm W-3. (11) 

Had we not made this subtraction correction, our limits 
on the cross section would be stronger by about an order 
of magnitude. 

C. INTEGRAL 

The space-borne INTEGRAL (INTErnational 
Gamma-Ray Astrophysics Laboratory) observatory [48| 
has searched for gamma-ray emission in the Milky 
Way over the energy range 20-8000 keV, using the 
SPectrometer on INTEGRAL (SPI). Teegarden and 
Watanabe [H| presented results of an INTEGRAL 
search for gamma-ray line emission from the Galactic 
Center region (we use their zero-intrinsic-width results, 
as appropriate for the low dark matter velocities of the 
halo). Other than the expected positron annihilation [5(1] 
and 26 Al decay [5l| signatures, no evidence of other line 
emission was found. 

To reduce backgrounds and improve the sensitivity of 
the line search, the measured intensity from large angular 
radii ( > 30°) was subtracted from that in the Galactic 
Center region ( < 13°), resulting in a 3.5-er constraint 
on the flux of very roughly < 10~ 4 photons cm~ 2 s" 1 
in the energy range 20-8000 keV. Our calculations must 
reflect this subtraction, which will somewhat weaken the 
sensitivity to the dark matter signal. A similar correction 
was used in Ref. [E^. We implement this as 

13° 

J An An = 2tt f (J(iP) - J{> 30°)) sin ipdip~2. 
Jo 

(12) 
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Due to the decreasing trend of the dark matter profile, 
the intensity outside the Galactic Center region will be 
largest at 30°, and accordingly we choose this value to 
be as conservative as possible (a larger subtraction leads 
to a weaker upper limit on (ajiv)). Had we not made 
this correction, our limits on the cross section would be 
stronger by about a factor of 2. 

D. Andromeda Halo Results 

The Andromeda galaxy (M31) has been observed by 
several gamma-ray experiments, all of which placed up- 
per limits on the flux. EGRET, CELESTE, and HEGRA 
all observed Andromeda, each encompassing a respec- 
tively smaller angular region of that extended object. As 
the results were reported as flux limits from specified an- 
gular regions, we compare to these using j7AfiAf2, which 
is an input to Eq. |J5]). 

EGRET viewed Andromeda with an angular radius of 
0.5° and set a 2-a upper limit on the gamma-ray flux of 
1.6 x 1(T 8 photons cm" 2 s" 1 from 0.1 GeV to 2 GeV, 
since no signal was seen (53|. For the angular region of 
this observation, the flux will be proportional to 

JaoAO = 2tt / J'(ip) sinipdip ^ 2 x 1CT 3 . (13) 
Jo 

CELESTE (CErenkov Low Energy Sampling and Tim- 
ing Experiment) is an atmospheric Cerenkov telescope in 
the French Pyrenees, which studies gamma rays with en- 
ergies greater than 50 GeV [54]. It viewed Andromeda 
in the en ergy range of 50-700 GeV, and again no signal 
was seen [551 ] . A 2-a upper limit on the energy- integrated 
flux from Andromeda was reported as < 10~ 10 photons 
cm~ 2 s — employing an angular radius of o bs = 0.29° 
yields J A nAfJ ~ 1 x 10~ 3 . 

HEGRA (High Energy Gamma Ray Astronomy exper- 
iment) was an atmospheric Cerenkov telescope, located 
in La Palma in the Canary Islands (56|. It took data in 
the range 0.5-10 TeV, with better energy resolution than 
that of CELESTE [57j. It used an even smaller angular 
radius of 6 ohs = 0.105°, which yields JaqAQ ~ 2 x 10~ 4 . 
HEGRA reported 99% C.L. upper limits for the gamma- 
ray line flux, and these can be used directly. 

E. Cosmic Diffuse Results 

INTEGRAL COMPTEL 59] and EGRET [U 

have all made measurements of the gamma-ray flux at 
high latitudes, and these can be used to set a limit on the 
cosmic dark matter annihilation signal. The INTEGRAL 
data used here were those collected in broad energy bins, 
much like those of COMPTEL and EGRET. The cos- 
mic gamma-ray background was also measured by the 
Gamma-Ray Spectrometer aboard the Solar Maximum 
Mission (SMM) [gj over the energy range 0.3 - 8 MeV, 
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m x [GeV] 

FIG. 2: The limits on the partial cross section, ((JAv)-y^ : 
derived from the various gamma-ray data. Our overall limit 
is shown as the dark shaded exclusion region. For comparison, 
the light-shaded region shows the corresponding limits for the 
NFW (rather than the Kravtsov) profile. 

for a field of view of 135° in the direction of the Sun [o2l ]. 
and we include this data. 

For the cosmic diffuse analysis, the framework detailed 
in Section IIII CI can be applied. Note that for simplic- 
ity we calculate only the true cosmic diffuse dark mat- 
ter signal, neglecting any Galactic contribution along the 
lines of sight. This contribution from the Galactic halo 
(which would add to the signal and thus make our lim- 
its stronger) is significant for NFW or steeper profiles 
and can even dominate over the true cosmic dark matter 
signal; see YHBA and Ref. 

V. DISCUSSION AND CONCLUSIONS 
A. Limits on the Cross Section to Gamma Rays 

In Fig. we combine all of the upper limits on the par- 
tial cross section to monoenergetic gamma rays, choosing 
the strongest limit for each value of the dark matter mass. 
The shaded exclusion region shows our combined bound. 
These searches for dark matter signals are limited by as- 
trophysical backgrounds, and the general trend of how 
the limits vary with mass follows from how these back- 
grounds vary with energy. We can estimate how the cross 
section limit should scale with mass, and how it should 
depend on the assumed spectrum of final-state gamma 
rays and the choice of density profile. 

Recall that we conservatively require the signal to be 
as large as the full measured background in an energy 
bin. The gamma-ray number flux of the signal inte- 
grated in a logarithmic energy bin A (mi?) scales as 
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Ed$/dE A(lnE) ~ (aAv)/m^, provided that the bin is 
wide enough to contain the full signal. The gamma- 
ray number flux of the background integrated in the 
same logarithmic energy bin scales as Ed&/dE A (In E) ~ 
E/E a A(\nE), for a background spectrum d<fr/dE ~ 
1/E a . For a narrow bin, the evaluation point is E ~ m x . 
We then expect the upper limit on the cross section to 
scale as (<JAv)u m it ~ w^ - " A(ln£'). For example, for the 
EGRET diffuse data, a is slightly greater than 2, and so 
the cross section limits in this energy range scale slightly 
less rapidly than as {<JA v )limit ~ m x- 

Most of these experiments had modest energy reso- 
lution. To be conservative, we assume an analysis bin 
with a logarithmic energy width of 0.4 in log 10 E (i.e., 
A(ln-E) ~ f) for the Galactic and cosmic diffuse analy- 
ses; this is at least as wide as the energy bins reported 
by the experiments. That is, even though we nominally 
assume two monoenergetic gamma rays at E 7 = m x , our 
results have not taken advantage of this fact. In effect, 
our results are what one would obtain for an annihila- 
tion gamma-ray spectrum as wide as 0.4 in log 10 E. The 
exception is the INTEGRAL line search, where the ex- 
cellent energy resolution is what leads to this limit being 
stronger than expected from the general trend in Fig. O 

Due to radiative corrections [27| or energy-loss pro- 
cesses [Hj], there should be some gamma rays near the 
endpoint, and our results can be scaled if the assumed 
branching ratio is less than the 100% used in Fig. [5J For 
example, for charged-particle final states, the branching 
ratio to internal bremsstrahlung gamma rays near the 
endpoint is Br(~f) ~ a ~ 10~ 2 . For neutral final states, 
there will typically be gamma rays (or neutrinos) near 
the endpoint. To be conservative about these details, we 
chose a nominal minimum branching ratio to gamma rays 
near the endpoint of 10~ 4 . 

How would our results change if we considered an even 
broader annihilation gamma-ray spectrum? We empha- 
size that the results shown in Fig. [5J which are based on 
direct numerical integration, are already valid for spec- 
tra as wide as our analysis bins. First, we should take 
into account the increase in the logarithmic bin width. 
Second, to be more precise, the evaluation point for the 
background spectrum should not be E = m x , but rather 
E = m x /a, with a > 1. This increases the estimate 
of the integrated background, and hence the cross sec- 
tion limit, by a factor ~ a" -1 . Thus, if we took the 
annihilation gamma-ray spectrum to be as much as one 
order of magnitude wide, then our limits in Fig. [2] would 
be weakened by at most a factor of several, depending 
on the background spectrum. (For the INTEGRAL line 
search, the correction would be much larger.) 

Given the large range on the axes in Fig. [51 and our 
intention to define approximate and conservative limits, 
this shows that our results are much more general than 
they first appear. Similarly, the results in BBM [7| and 
YHBA [8| do not have a strong dependence on assumed 
annihilation neutrino spectrum. 

How sensitive are our limits to the choice of density 



10 




rn, [GeV] 

FIG. 3: The gamma-ray and neutrino limits on the total an- 
nihilation cross section, selecting Br^^j) = 10 -4 as a conser- 
vative value. The unitarity and KKT bounds are also shown. 
The overall bound on the total cross section at a given mass 
is determined by the strongest of the various upper limits. 

profile? As noted, we chose the rather shallow Kravtsov 
profile to be conservative. If we were to adopt an NFW 
profile, which increases much more rapidly toward the 
Galactic Center (scaling with radius as r" 1 rather than 
r~°' 4 ) the annihilation rates would be larger and the cross 
section limits correspondingly stronger. In Fig. [2j we 
show how our results would change if we had used an 
NFW profile instead of the Kravtsov profile. At most 
energies, the changes are modest, and illustrative of the 
potential uncertainties. The only significant change to 
the combined gamma-ray limit is for the H.E.S.S. Galac- 
tic Center Ridge case, which is based on small angular 
radii. In the NFW case, the steeper profile gives an over- 
all larger intensity and a smaller signal cancelation when 
the background is subtracted. A fuller discussion of how 
the annihilation signals depend on the choice of dark mat- 
ter density profile is given in YHBA. 

B. Limits on the Total Cross Section 

Unsurprisingly, the cross section bounds derived under 
the assumption of monoenergetic gamma rays are sub- 
stantially stronger than those defined similarly for final- 
state neutrinos in BBM [7] and YHBA 8] . (At the high- 
est masses, near 10 4 GeV, this is no longer true, first be- 
cause of how the numerical limits work out, and then be- 
cause we do not presently have good gamma-ray data or 
limits at higher energies; we expect that dedicated analy- 
ses by H.E.S.S. and other experiments will soon improve 
this.) Indeed, this was an assertion in those two works 
that we have now justified in more detail than before. 
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It is unrealistic to have Br (77) = 100%, of course, if 
one is trying to set a limit on the total cross section. If 
Br(77) is known, then a limit on the total cross section 
can be determined by dividing the limit on the partial 
cross section to that final state by the branching ratio: 



(o-av) 



total 



BK77) 



(14) 



In typical models, this branching ratio is typically 10 -3 or 
smaller [lj, 2, 3]. To be conservative, we must just choose 
a value such that it is implausible that the true branching 
ratio could be smaller. We therefore assume Br (77) = 
10~ 4 , but this choice could be debated. As noted, our 
analysis uses wide logarithmic energy bins, and so, at 
the very least, would capture the gamma rays near the 
endpoint due to internal bremsstrahlung from charged 
particles [27] . (Similarly, as a general point, limits on 
the total cross section defined by assuming only W + W~ 
final states [HJ would have to be corrected by dividing 
by Br(W+W-).) 

Figure [3] summarizes various limits on the total cross 
section, including the one just described, the unitarity 
bound mentioned earlier, and the neutrino bound from 
YHBA (based on the Milky Way signal and the Kravtsov 
profile) . The standard cross section for a thermal relic is 
also shown. Note that our limits bound (<jav) directly, 
independent of whether a a is s-wave or p-wave domi- 
nated. These results, combined with those in our Fig. [21 
strongly constrain the possibilities for large dark matter 
annihilation signals, e.g., as assumed in Ref. [64| . 

When shown in this way, it becomes clear how surpris- 
ingly strong the neutrino bound on the total cross sec- 
tion is, as it is comparable to the bound obtained using 
the gamma-ray flux limits and a reasonable assumption 
about the minimum branching ratio to gamma rays. It 
is very important to emphasize that while the gamma- 
ray bound on the partial cross section had to be divided 
by a realistic Br (77), this is not the case for the neu- 
trino bound, as explained above. If we assume only SM 
final states, then all final states besides neutrinos lead to 
appreciable fluxes of gamma rays, and hence are more 
strongly excluded. Of course, the gamma-ray and neu- 
trino cross section limits can both be weakened by as- 
suming an appreciable branching ratio to new and truly 
sterile particles. 



C. Conclusions and Prospects 

Using gamma-ray data from a variety of experiments, 
we have calculated upper limits on the dark matter an- 



nihilation cross section to gamma rays over a wide range 
of masses. These limits are conservatively defined, in 
terms of our analysis criteria, our assumptions about the 
uncertain dark matter density profiles, and the gamma- 
ray spectrum. While our results were nominally defined 
for monoenergetic gamma rays with B 7 = m x , we have 
shown that all of our results except the INTEGRAL line 
flux limit are only weakly dependent on this assumption. 
The limits obtained for more general gamma-ray spectra 
would only be somewhat less stringent. 

There are good prospects for improved sensitivity with 
present and upco ming gamma-ray experiments, particu- 
larly GLAST [gE II Ed and the TeV ACT detectors. 
More detailed searches and analyses by the experimental 
collaborations themselves should also lead to improve- 
ments, which we encourage. These searches for dark 
matter signals are already background-limited, which will 
limit the possible improvements. GLAST and other ex- 
periments should be able to make reductions in the back- 
grounds by taking advantage of better energy and angu- 
lar resolution, and by reducing the residual diffuse emis- 
sion by subtracting astrophysical components and resolv- 
ing individual sources. The high statistics expected for 
GLAST and other experiments should also make it possi- 
ble to define detection criteria in terms of the uncertainty 
on the background, instead of the whole measured back- 
ground. 

Using a conservative choice on the branching ratio to 
gamma rays, namely Br(-yj) ~ 10 , we defined an up- 
per limit on the total dark matter annihilation cross sec- 
tion by dividing our limits on the partial cross section 
to gamma rays by this branching ratio. At intermedi- 
ate energies, the upper limit on the total cross section 
defined this way is comparable to previous upper limits 
defined using neutrinos [3, El, . The combined limit is 
considerably stronger than the unitarity bound [24], I25I ] , 
or the cross section of Ref. [|| , which would lead to sub- 
stantial modifications of dark matter halos. For the rela- 
tively large cross sections considered here, the dark mat- 
ter could not be a thermal relic; additional work is needed 
to push the sensitivity of these and other techniques down 
to the expected cross section scale for thermal relics. 
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